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Magnetostatic  Surface  Wave  Transducers 


Introduction 

The  purpose  of  this  report  is  to  summarize  the  results  of  the 
work  on  magnetostatic  surface  wave  transducers  under  contract  number 
F  19628 -80-C- 0029  from  the  U.S.  Air  Force  ESD  RADC  EEA  at  Hanscom 
AFB,  Ma. 

To  be  presented  are  the  results  for  the  dispersion  relation, 
radiation  resistance,  radiation  reactance  and  insertion  loss  for 
magnetostatic  surface  wave  transducers  which  may  include  a  gap  and 
apodization.  Independent  conductors  as  well  as  normal  modes  are  con¬ 
sidered.  Also  presented  are  the  results  for  the  dispersion  relation 
for  surface  waves  for  a  variety  of  alignments  of  the  externally 
applied  magnetic  biasing  field. 

Computerized  results  including  computerized  graphs  of  the  results 
are  presented  here.  Comparisons  are  made  with  results  obtained  from 
the  analysis  of  the  microstrip  model  which  is  also  here  presented. 
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Basic  Theoi 


The  basic  theory  leading  to  the  dispersion  relation,  magneto¬ 
static  wave  power,  radiation  resistance,  radiation  reactance  and 
insertion  loss  for  surface  waves  when  the  applied  magnetic  biasing 
field  is  in  the  direction  of  the  Z  axis  (see  Figure  1)  has  been  pre 
viously  detailed  ^  ^  ^  ^  .  This  theory  will  here  be 

outlined. 

We  start  with  Maxwells  equations 


m.§ 


VxE  = 


V-B  -  o 


V-D  =  0 


and  the  constitutive  relations  in  each  of  the  four  regions 
B  -  uQ  (H  +  H) 

(2 

D  =  eE 

where  M  is  taken  as  zero  in  all  regions  except  the  YIG  region.  We 
utilize  the  gyromagnetic  relation  in  the  YIG  region 


=  -Y  M  x  H 


and  retain  only  first  order  terms. 

We  assume  the  time  dependence  of  all  physical  quantities  to  be 
e^wt.  We  also  take  the  magnetostatic  approximation 


H  =  E  =  E  =o 


Si 


and  no  variation  of  any  physical  quantity  in  the  z  direction. 
In  particular,  we  obtain 

3Ez 

“17  =  '  J  i0  Bx 

3E 

"if  =  +  j  w  Bv 


(5) 


and 


E  =  u  H 

X  OX 


B  =  u  H 
y  o  y 


(6) 


in  all  regions  except  the  YIG,  while 


'JU12 

y/  Vjy21  y22, 
in  the  YIG  region  where 

M  =  V 


11  22 


=  1  - 


H 


__ - £ 

n  - 


n 


u  _  y  .  _ 

21 '  12 " 

a  =  JLl 

4ttM 

o 


<Vi  = 


H  OT 


(7) 


(8) 


Y  =  2.8  mhz/oe  ;  4ttM0  =  1750  oe 
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f  =  co/2tt 


Solutions  are  sought  which  satisfy  continuity  conditions  for  Hx 
and  By  at  each  region  junction  and  satisfy  By=o  at  the  ground  planes. 
At  y=g  the  condition  to  be  satisfied  is  that  is  discontinuous  by 
the  surface  current  density  J(x). 

We  thus  assume  a  solution  form  of  a  potential  function 

4>  =  F(y)  e^'K*)  (9) 

where 

Hx  ■  S  *  >V  -  $  (10) 

In  the  non  YIG  regions  we  find  the  form  of  F(y)  to  be 

F(y)  =  Ai  e!k'y  +  Bi  e'^y  i-1,3,4  (11) 

while,  in  the  YIG  region 

F(y)  =  A2  e0lk!y  +  B2  e'^'y  (12) 

where 

e2  =  yll/y22  (13) 

so  that  the  basic  equations  (1)  -  (8)  are  satisfied.  One  can  see  that 
these  solutions  consist  of  waves  propogating  in  the  X  direction.  We 
carry  8  along  in  the  analysis  even  though  its  value  is  unity  by  (13) 
and  (8)  because  of  comparisons  to  be  made  later  with  the  analysis  for 
a  general  direction  of  the  applied  biasing  field. 

The  attenpt  to  satisfy  the  continuity  and  boundary  conditions 
results  first  in  the  requirement  to  solve^ 
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Ft(K)  =  o 


(14) 


where 


(K)  =  ®-coth I ^ I ti ' 1 )  f Cl+a^)e'2e'k'd+Cl-an)T]e''K'g  (cothlKlt1+1) 

1  2  2  ~ 

[(l-a2)e‘26^Kld+(l+a1)Tle^Kfg  (15) 


and 


al  =  y22  6  +  V*  u12 


a2  y22  6 


M 


K  y12 


T  =  (a2+tanh  j  K | £) 
"(u^-tanh  ]K].5.) 


(16) 


Equation  (14),  a  transcendental  equation  for  K  as  a  function  of  f, 
is  the  dispersion  relation.  Numerical  techniques  are  required  for  its 
solution.  Two  solution  curves  of  K  vs.  f  result;  in  one  solution  K  is 
always  positive  and  in  the  other  solution  K  is  always  negative.  This 
results  in  two  solution  waves  which  are  in  opposite  directions.  Denoting 

^  =  s  (1?) 


and  the  solution  values  of  K  by  Kg,  S=-l,  1,  we  have  that  the  two  dis¬ 
persion  relation  curves  are  obtained  by  solving  (14)  with  (16)  for 
S=-l  and  S=l. 

Equation  (14)  can  also  be  written  as^ 
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(18) 


e-2|K|T 


(l-a2)e'2BlK'd+(l+a1)T 

(l+a2)e'23,Kld+(l-a1)T 


where 


x  =  +  g 


(19) 


which  shows  that  the  effects  of  material  thickness  t^  and  g  enter  the 
dispersion  relation  only  in  combination. 

(131 

Another  useful  way  of  writing  the  dispersion  relation  isL  J 


e~26 1 K | d  (ctj+tanhjKjx)  (a2+tanh|K|£) 

~(a 2 - 1 anh [ K | t )  (ot^- tanh[K 1 1) 


(20) 


'Ihe  bandwidth  of  frequencies  for  which  the  solution  of  (14)  can 
be  obtained  is  given  by^ 


yVho(Hq+4ttMo)  <  f  <y  (H0+2ttMo) 


(21) 


Having  the  dispersion  relation  curves  we  can  find  the  group 

delay 


V 


g 


3(a) 

3K 


(22) 


for  each  of  the  two  solution  curves. 

After  equation  (14)  has  been  solved  we  can  find  all  quantities  of 

|2i 

physical  interest  for  each  of  the  two  solutions  of  (14)  .  The  mag - 

r 21  (31 

netostatic  wave  power  is  then  obtained  from  1  1 ’ 1  J 

p(s>  ,  I  f  ^  dy  S--1.1  (23) 

-oU 

where  is  related  to  H^,  and  Hx  in  the  YIG  region,  by  equations  (5) 
with  (6)  or  (7) . 
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The  expression  for  power  is  found  to  be 


[2] 


>(s)  =  -su)yQ  A  g2 
- s  s 

2  r 


s=-l,l 


(24) 


where 


G  = 
s 


e'^sl^oyl 


la|-FT(K)l 


K=K 


A  =  Cyi) 

S  - 55- 


s«-l,l 


sinh  2 1  Ks  J  £  |Ks|£\  (Ug  elKs^Vs  e'lKs'8) 

cosh"  I  Ks  |  £  y  4  ‘  -I-)*  4~sinh^ 


(25) 


[sinh  2 S Ks ! t1 


|Ks!ti 


ru2 


4  (e2,Kslg-l)-  -J  (e*2iKs'g-2)-2  Ug  Vs  |K Jg 


a(s)  t2  (e26'Ksid-l)-a^-)(e'2BlKs'd-l)-2  B2  |Kc|Tcd 


s  s  22 


(26) 


S=-l,l 


Us  =  +  (l+ajs^)Ts  eBlKs'd 

Vg  =  (l+a^)e'6'Ks'd  +  (l-ajs))Ts  e6'Ks'd 

S-1,1  (27) 

“{S)  ■  “l«y  .  «?’  -  “2(V  -  Ts  ’  T(Ks) 


For  independent  conductors ^ * ^2^ 


N 


J,  (K  )  =  Z  sine  aiKs  ttYE-  e^  Ks  Pi  1 
1  S  i=l  2tt  11 


-  7  - 


s=-l,l  (28) 


For  the  non-apodized  independent  conductor  case,  (28)  can  be  written 
as  ,  with  I  =1, 


a  K  l-nN  ejKspN 
Jl(Ks)  =  IQ  sine  ^s  1  n-  e  s 


1-n  e^Ksp 


(29) 


For  a  truncated  array  of  normal  modes  we  have  for  the  fundamental 
mode  (n=l) 


N 

J,  (K  )  =  E  sine 
1  s  i=l 


2  a. 


Pi(3-n) 


sine 


K  p. 
s1  1 

2tt 


3+n 

~r 


niV^"  e-i  Vi1 


(30) 


s— 1,1 

where  a^,  p^  i=l,2,...,N  are  the  conducting  strip  lengths, 
conducting  strip  widths  and  center  to  center  spacings,  respectively, 
to  account  for  apodization.  N  is  the  number  of  conducting  strips  and 
n=-l  for  a  meander  line  and  n=l  for  a  grating. 

The  definition 


sine  X  =  sin  ttX 
ttX 


(31) 


is  employed  in  the  above. 

In  the  free  space  case,  t^=°°  and  £=°°,  the  dispersion  relation  is 


e  -  2  6 1 K  |  cl  =(oti+1)(a2+li 
(a2-l)  (a1-l.) 


(32) 


from  (20).  The  expression  for  Ag  in  the  power  is 


V(T 


1s>1)  +  ^  T2  (e2BlKs'd-l)-  *2^  (e~2BlKsld-l)-282Tsvi22d| 


il  V2 

K  I  +_s_ 

J  8 


(33) 


8 


The  radiation  resistance  is  then  obtained  from 


R(s)  =  4  jP(s)l  s— 1,1  (34) 

(l-n)+(l+n)N2 

In  the  above  it  is  typical  that  the  amplitudes  for  the  wave  corres¬ 
ponding  to  s=-l  is  greater  than  the  amplitudes  for  the  wave  corresponding 
to  s=+l.  Thus  the  s=-l  wave  is  the  stronger  of  the  two  and  is  denoted 
by  the  +  wave  and  the  s=l  wave  is  denoted  by  the  -  wave. 

The  total  radiation  resistance  is  then 

R  =  R+  +  R~  (35) 

m  v  ' 

The  radiation  reactance  contributes  meaningfully  for  surface  waves. 
It  is  to  be  obtained  from 

CO 

xm  mh  f  Rh,(£,)  d£'  (36) 

m  71 J  F~T 

-oo 

Although  this  integral  contains  infinite  limits  and  an  apparent 
singularity  it  can  be  computed  accurately  by  numerical  techniques  ^ ^ 
We  thus  find  the  radiation  reactance  from  the  numerical  values  of  the 
previously  obtained  radiation  resistance. 

The  combination  of  radiation  resistance  and  radiation  reactance 
then  results  in  the  complex  radiation  impedance. 

We  can  now  find  insertion  loss  from  the  radiation  resistance, 
radiation  reactance  and  by  including  source  resistance,  conduction  loss, 
matching  reactance  and  propagation  loss.  We  obtain  ^ 

IL(s)=20  log, n  4r(S:>  Rg  _  76.4  AH  Ar  s=-l,l  (37) 

~ssm 
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where  R  is  the  source  resistance,  R,  is  the  conduction  loss  and  X, 
g  L  L 

is  a  series  matching  reactance.  AH  is  a  linewidth  representing  material 
loss  and  Ar  is  a  propagation  distance. 

This  completes  the  basic  theory  for  magnetostatic  surface  wave 
transducers . 


Microstrip  Model 

In  the  microstrip  model  insertion  loss  is  calculated 

from  input  resistance  and  inputt  reactance  of  a  lossy  shorted  section 
of  microstrip  line  and  microstrip  propagation  constants.  Apodization 
is  not  taken  into  account  in  this  theory.  The  conducting  strip  dimen¬ 
sions  fcp  a  and  p  are  thus  the  same  for  all  strips. 

First  consider  one  conducting  strip.  Here 


J(Kg)  =  sine  e~jKsP 


(38) 


from  the  independent  conductor  model  ,  (.28) . 

For  N  strips  we  then  have  ,  multiplying  by  an  array  factor, 

NKgp  . 2 


for  n=l, and 


,(s) 


4|P(S)| 

=  u^nT+C 


sin 


1+n) 


T 


K  p 
sm  s 


s—1,1  (39) 


R(s)  =  4jPCs)( 

(l-n)+(l+n) 


~T 


NKgP 

sin  — 2~ 

COS  -J~ 


s— 1,1  (40) 


for  n®  -1  and  N  even  where,  is  computed  for  N=1  in  (39)  and  (40), 

From  the  above  and  (35)  and  (36)  we  have  R+,  R  ,  and  X^. 
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To  obtain  insertion  loss  we  first  define 


7^- 

1/2 


h/2 

h/2 


s— 1,1 


(41) 


Given  characteristic  impedance  Z propagation  constant  Bc, 
conduction  loss  constant  a  ^  and  conductivity  o,  we  have,  for  N=1 


“R 


a 

c 


(42) 


where  and  ac  represent  radiation  attenuation  loss  and  conduction 
attenuation  loss,  respectively. 

For  one  or  more  conducting  strips,  N>1,  we  now  have 


c 

“R 


Icl/N 


*  YN 


6  =  6, 


a_  =  a 


nBl 


(43) 


and 
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z_  =  z 


cl 


“r  ■  VN 


n=-l  and  N  even  (44) 


8  = 


8X  N 


a 


a  N 
c 


Total  attenuation  loss  is  then 


0  3  *R**c 


(45) 


Then  compute 


R. 


Z„  tanh  2  a  JL 
c  l 


in  1+cos  2  tftj  sech  2  a  ^ 


X, 


Zc  sin  2  15  ^  sech  2  a  ^ 


in  1+cos  2  sech  2  a 


in  V  i 


2  +  X-2 
m  in 


(46) 


where  R-  ,  X.  and  Z.  are  input  resistance,  input  reactance  and  the 
in  in  in  r 

magnitude  of  the  input  impedance,  respectively.  With 

s— 1,1 


R^S)  = 

i,m 


R. 

in  — - 


VRm  '  Zc 


(47) 


R.  =  R.  +  R.  m 
i,m  i,m  i,m 


insertion  loss  can  be  expressed  as 


IL<S>  =  20  log10  4  Rj  R^> 

- 7 


76.4  AH  Ar 
3w/  3K 


s— 1,1  (48) 


(R.+R.  )  +  X-  ‘ 

v  i  in'  m 
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We  here  indicate  the  computation  of  the  magnetic  wave  power  P  when 
all  physical  quantities  are  determined  by  combining  the  two  solutions 
present.  We  shall  only  consider  the  free  space  case  (il=tj=«0  and  the 
case  of  no  gap  present  (g=o) .  Since  6=1  in  the  basic  theory,  we  will 
eliminate  it  from  the  equations. 

We  write 


(  s') 

°1  ‘  y22  *  S  v12 
“2  '  *  v22  '  S  u12 


s=-l,l 


(51) 
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and  define 


The  magnetic  wave  power  defined  as 


is  now 


op 

P  -  \  f  Ez  Hy  dy  (55) 

-CO 

r  -d  o  » 

■  7  / Ez1  Ryl  **  *  f  Ez2  Hy2  *  *  /  Ez4  V  W 

L  -oo  -d  n 


Considering  both  solutions  we  have 
d 


P4[/  (E2Ci1)*Ezl))(H^i1)*H^l))dy*/  (^21)+Ez12))a)y21)+Hf2):idr 

/  <E£  t55) 

o  -* 

Employing  (5)  with  (6)  or  (7)  together  with  (9) -(12)  indicates  the  fonn 


of  (55)  is 


jj  |k^)e’jK-lX  +  E1(y)e';iKlxJ[Hy_i(y)e^K-lx  +  H^V]  dy 

\  -00 

+  /  [E.2(y)e"jK-lx  +  E2(y)e'jKlx][Hy_2(y)e^K-lx  +  H^Cy^V]  dy 
-d 

*  /  [E.4(y)e-JK-lX  *  E4(y)e-)ltl!‘][Hy.4<y)e)K-lX  *  ^ 
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and,  upon  insertion  of  the  appropriate  solution  functions  in  (55)  and 
the  performance  of  the  indicated  integrations  we  obtain  the  result 


P  =  ^o  ^Mg  +  Mg  cos  (Kj-K_1)x  +  j  My  sin  (K1-K_1)x| 


(57) 


where 


1  s  Gs2 

V  s  — T 

y  s=-l  2  K 


(1+a/  +  (e2!Kstd-l)-a^  (e'2!^-!) 


- 2 1 Ks | as  (ajs)  +  a^^d  +  (as  e^s^d-a|s^  e^s^d) 


(58) 


Mg  -  TO  (Ua.pa.ap^^a^e-tV^l^-n-a^a.^e^-l^- 

'^i  K_i  |  (K^-K _^) 


♦  (e(Kl'K-l)d-l)C(tj1)K14pK  j)-a  jSj  (e‘ (ltl'K-l)d-l)  x) 

- - Tifprp -  (59) 

*  dCj^.p  (a_j  e-|t-ld)(h  S?  eh^eh*)] 

M,  aK-ld-^«e-K-ld)(a1  e‘ V-a^eV) 


TkT^Tc 


h(e-¥.l)d.l).,«  a.1(a«l*i:-l)d-l)] 
(e^i'K-i^-i)^"1)  K1-^)K_1)-a_1a1(e"(Kl'K-i:)d-l)(a2('1^K1-a^1)K  x) 

CTO 


(60) 


Note  that  (57)  is  of  the  form 

P  =  PR+jPl 


(61) 
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s 


where 


-  [Mg  +  Mg  cos  O^-K.^x 


Pj  =  ~2 —  sin  (K^-K _^)x 

The  results  obtained  reduce  to  that  obtained  earlier  in  (24) , (25)  and 
(33)  when  the  two  solutions  are  considered  separately  and  (50)  and 
the  dispersion  relation  (31)  is  utilized.  In  this  case  the  only  terms 
present  in  (57)  are  the  two  terms  in  Mg,  one  for  each  value  of  S,  from 
(58). 

In  general,  the  complex  impedance  is  taken  as 

Z  *  - L_P -  (63) 

(l-n)+(l+n)Nr 

similar  to  (34)  and  has  real  and  imaginary  parts.  The  spatial  average 
of  this  generalized  impedance  gives  the  radiation  resistance,  while  the 
spatially  dependent  part  gives  rise  to  resistance  and  reactance  terms 
related  to  the  width  of  the  transducer  in  the  x  direction.  These  terms 
are  assumed  to  be  of  second  order  and  have  not  been  incorporated  into 
the  present  model. 

Generalized  Dispersion  Relation 


In  this  section  we  obtain  the  dispersion  relation  for  surface  waves 
when  the  biasing  field  is  not  restricted  to  be  parallel  to  the  Z  axis 
(see  Figure  1) . 

In  the  YIG  region  the  components  of  the  permeability  tensor  (7)  are 
now  given  by 
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yll = 
y22  = 
■jy12  = 
jy21  = 


1  + 


1  + 


y^H  (4itM  )  (sin20  sin2<J>  +  cos20) 

2  u2  c2 
Y  H  -£ 

Y^  (4ttM  )  sin20 
o  o 

— 2  772  72 

Y  H  -f 


j  y(4ttM  )sin6  (£  sin0  +  jYHQ  cos4>  cosB  ) 

Y  H  -£ 

-j  Y(4irMo)sin0  (f  sin0  -  jYHQ  cos<f>  cos0) 

2  u2  ,2 

Y  H  -f 


(64) 


These  relations  reduce  to  those  given  by  (8)  for  the  case  of  the  biasing 
field  lying  along  the  z  axis;  0=90°  and  $=90°. 

The  satisfaction  of  (l)-(3)  yields  solutions  as  in  (8),  the  ex¬ 
pression  in  the  YIG  region  being  modified  to 

F(y)  =  e';jKby(A2  e6!Kly+B2  e'e'Kly)  (65) 


instead  of  (12). 


Here 


S  = 


(y2l‘y12)  +  4yll  y22 


>  o 


4p 


22 


(66) 


instead  of  (13),  and 
b 


"j (y21  ~  y12j 
2  y22 


(67) 


We  note  that  b  is  real  and  that  the  term  containing  b  indicates  that 
there  is  an  additional  propagation  component  in  the  y  direction.  For 
the  standard  surface  wave  case  of  0=90°  and  <j>=90°  we  note  that  (65) 
and  (66)  reduce  to  (12)  and  (13)  with  b=o. 

The  dispersion  relation  is  obtained  by  satisfying  the  boundary 
and  continuity  conditions  as  before.  With 
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(68) 


“l  =  0  ^22  ‘  j  ^  U21  +  b  V225 

a2  =  6  U22  +  3  ¥■  y21  +  b  W22} 

the  dispersion  relation,  for  the  case  of  g=o,  is 

e2g|K|d  _  jK|t^) (a^-tahh  [ K [ £.) 

(a2+tanh  |K|0  (a^+tanh  ( K( t^) 

Again,  for  the  case  of  0=90°  and  <£=90°  equations  (68)  coincide  with 
(16)  and  (69)  is  then  the  same  as  (20) . 

Thus  (69)  with  (68),  (67),  (66)  and  (64)  give  the  dispersion  rela¬ 
tion  for  surface  waves  with  the  orientation  of  the  biasing  field  kept 
arbitrary. 
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COMPUTER  PROGRAMS 


A.  Basic  Theory 

A  computer  program  which  incorporates  the  results  of  the  basic  theory 
has  been  made  operational  on  the  CDC  6600  at  Hanscom  AFB,  Ma.  The  program 
produces  plots  of  the  various  physical  quantities  as  functions  of  frequency. 
There  are  plots  of  wave  number,  group  delay,  radiation  resistance  and 
insertion  loss  for  each  of  the  two  solution  waves.  There  are  also  plots 
of  the  normalized  dispersion  for  the  +  wave,  total  radiation  resistance 
and  the  corresponding  total  radiation  reactance.  The  program  also  pro¬ 
vides  for  print  out  of  these  quantities. 

The  program  is  designed  for  flexibility  in  that  independent  conductors 
as  well  as  a  truncated  infinite  array  of  normal  modes  can  be  accomodated. 

The  case  of  uniform  conducting  strips  can  be  handled  as  well  as  apodization 
in  strip  length,  strip  width  and/or  center  to  center  spacing.  In  addition, 
the  program  automatically  computes  the  relevant  frequency  range  by  utili¬ 
zing  (21). 

There  now  follows  a  detailed  description  of  the  input  cards  to  the 
program  which  shows  how  to  use  the  features  described  above. 

Card  1  -  HQ,  tj,  d,  g,  N,  n 

These  seven  quantities  are  here  supplied,  separated  by  commas.  All  lengths 
are  in  meters.  Colunns  1-72  may  be  used. 

Card  2  -  first  ,  Afc^,  £,  option 
Card  3  -  first  a,  Aa,  a  option 
Card  4  -  first  p,  Ap,  p  option 
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Each  of  the  above  three  cards,  applying  to  Jl^,  a  and  p,  respectively, 
contain  three  items,  separated  by  commas.  The  first  item  is  the  di¬ 
mension  of  the  first  strip.  If  the  third  item  (option)  is  o  then  the 
dimensions  of  the  (N-l)  strips  following  the  first  are  successively 
incremented  by  the  increment  (A)  of  the  second  item.  If  the  third 
item  is  1  then  the  dimensions  of  the  first  —  strips  following  the 
first  strip  are  successively  incremented  by  the  increment  of  the  second 
item  and  the  dimensions  of  the  next  strips  are  successively  de¬ 

cremented  by  the  same  value.  N  must  be  odd  when  the  option 
is  1.  Note  that  one  may  use  a  negative  number  for  the  increment  of 
item  2.  Also  note  that  the  quantities  5.^,  a  and  p  are  handled  entirely 
independent  of  each  other.  If  an  increment  value  is  o  then  there  is 
no  apodization  in  the  corresponding  quantity  and  the  option  is  immaterial. 
Card  S  -  AH,  Ar 

These  two  quantities,  separated  by  a  cornua,  are  here  supplied. 

Card  6  -  heading 

Oily  columns  1-20  are  used  for  this  card.  For  the  normal  modes  case, 
the  first  ten  columns  on  this  card  should  contain  NORM_MODE_.  For  in¬ 
dependent  conductors ,  the  first  ten  columns  should  contain  IND  _COND_. 

This  card  serves  as  the  top  heading  line  on  the  plots  as  well  as  to 
signify  the  program  whether  the  case  is  one  of  normal  modes  or  uniform 
conductors . 

Card  7  -  heading 
Card  8  -  heading 
Card  9  -  heading 
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These  three  additional  heading  cards  are  required  and  will  appear  in 
order  on  the  computer  plots  under  the  heading  of  card  6.  Columns  1-70 
may  be  used. 

There  follows  a  listing  of  the  entire  computer  program  as  it  is 
used  on  the  CDC  6600  at  Hanscom  AFB,  Ma.  Omitted  are  the  required 
control  cards  consisting  of  the  standard  job  card,  Fortran  compile  and 
execute  cards  and  the  standard  control  cards  for  plots. 


PROGRAM  ROOT  (INPUT, CLTPUT) 

DIMENSION  PROFIO(3) 

niHRNSION  F  ( 1  2  0  d )  ,FM(12G0)  ,FP(i2QG>  ,  CAP (  120  0)  ,CAM(12QQ),VGM(12dG>, 
VPP{1?OO),PM(12O0>,RP(12OO>,RM(12OO),RT(12OO),PX(12OO),SERP(12OO), 
XSFRM(1200) ,VNN  (1200) ,VGP<  12  301 
DIHEN5I ON  HPAf ( 2)  ,HFAC1 (7) , HE AD2  (?)  ,HEAC7 (7) 

PTMfc  NSI ON  FN(50), VM  (  5  0  > 

COMMON  FL  ,  At  1  ,AL2  ,8,C,T1  ,G,S,  ET  A  ,EN  ,P,AY,  A,  EL  1  (40), PE  (40)  ,  A  A  (40) 
X.LMOOE 

R  FA  n  %h,T1,C,G,EL,FN,ETA 
Rr  A  C  ♦, EL8EGN,ELDEL,EL0PT 
PFAD  ♦, APPGIN,ADEL,ACFT 
PFAO  +, PRFGIN,PDEl,POPT 
PtA  P  «,CIELH,CIST 
PFAD  10  2  »  Hf  AO 
REAP  10  0 , HEAD  1 
°F  A  P  101, HFA  C  2 
PFAD  101,H£AD? 

L  MQDF  =  1 

TP  ( HE  A  0 ( 1  )  .FO.  "NORM  MOOE  " I  LMOCE=2 
N=E  N 
PL=  30. 

PL=T<=. 

PO  41  T=1,N 

FL1(I)=ELPEGN*(I-1)*EL0EL 
AAi T) =A ^FGTN>  (I-l)*AOEL 

41  PF(I  )  =P=!FGIN  + <I-l)*POEL 
NFL  = ( N  + 1  )  /? 

IF  (FLOPT  .FO.  0.)  GO  TO  42 
00  4?  I=NFl,N 

43  Fl1(I)=FL1 <NFL)-(I-NEI)*EIDEL 

42  TF  (  AO PT  .EQ.  c.)  GC  TO  44 
00  45  T=NF_L,N 

45  AA(I)=AA(NFL>-<I-NEL>+ACEl 

44  TF  (PCPT  . FO.  0.)  GO  TO  4b 
OP  47  I=NCL.N 

47  PF(T)=Pf (NFL)-(I-NEL)*POEL 
4ft  CONTTNUr 

FLO=2.b*SOPT (  H* (H+175U.) ) 

FHI  =  ?,P‘  (  H  +“  7  c  •  ) 

FOFL  =  1 . 

F°EG  =  AINT  (FI  0+1. 

NF=TNT(FHI)-INT(FL0)-1 
no  uo  i=i ,  nr 

40  F (I ) =F^FG+  (T-  1  >* FPEL 

IF  (Ml)  .LT.  FLO  PRINT  * , "PR EQU ENCI ES  TOO  LCW" 

IF  (F(MF)  .GT.  FHI>  PRINT  *  ,  **  F  RE  CU  E  N  C  IE  S  TOO  HIGH" 

RRTNT  60 
FACT=1 . 

TF  (FT A  .GT.  -2.)  GO  TO  4 
+'  T  A  =  - 1 . 

FACT=(? . /ENI * -2 

P=1NT  * PI  GRATING  CASE" 

4  F  ON  T I  Nil  F 

PRINT  6 1,H,T1  ,n,G»EL  ,EN, £TA,NF 

PRINT  ncL  T  A  H  =  "  ,  0  £  L  H ,  DISTANCE  =",DIST 

PRINT  *,  ••  c  LOSS  IS  "  ,  RL 
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POINT  81,  «T,?uaM*t,H) 

PRINT  82,  f  I*  «0  (I) 

PRINT  63,  <I,PF  U)  ,T=1,N) 

PRINT  60 

TP  <FTft  ,rO.  -1.1  ELL  =.87*.4E-8*EN 
IF  (FT  A  .fp.  l.)  ELL".23*.4E-8/£N 
ELL=0. 

OAT A  PPOGIO/PHSf THAPES, 4H372* ,10HJ. WEINBERG/ 
CALL  FLTIO?(PeOGID»2QO.,i2.,i»l 
P=0  . 

As-I  • 

AY^G. 

°I=  3. 14  15‘>?65  4 
°G=50. 

ATT  -  .5*  AT*  ( ( 1  .-FTAH-ENMl.  +  ETAll 

1)0=4  ,*P  T*  1 
K -1 
1=1 
J=i 

OMH  =  H/l  7«5''. 

5C  FF-F(K) 

DM=FF/(2. P  *1 7F0  -  > 

Ult =  i.-OMW/f 
U22=U1 1 

U12=C*/ (OK*‘?-C!NH**2> 

P=S0PT(U11/U2?> 

L  =  1 

30  TP  (l  .  EO.  2»  GC  T  ? 

!  Ssl. 

GO  TO  3 

2  S=-i. 

3  C ONT  T NUr 

At  1  =  U2?*B»S*U12 
Al?=U22*R-S*Ul? 

TP  (J  .  GT,  11  GO  TO  53 

C  AO  =  ,5*OQ°T<U<-2/Ull)  *ALCG(1  .  *4  .  *  SQR  T  (L 1 1*  l>2  2)  / 

t(U12**2-(SO°T  <U1  l*(/2 Z)  *t,  1**21 ) 

CAO=CAO/D 
53  CONTINUE 

ip  n  .  fo.  i>  r,r  to  51 

IP  <J  .FO.  11  GO  TO  51 
IP(l  .EO.  11  GACsCAPU-1) 

IP (L  .€0.  21  C  AO=CAM <  J-l  I 
51  Hsl 

5  PPL  =  .02*  CAP 
C  AOPsCAO+Cft 
C  A OM=C A  O-OPl 
CAPO=CAO*  D 
C  ADD  =  ABS ( C  AOO  ) 

CAOGrAPS(CAO‘G) 

IF1CACD  .GT,  650.1  GC  TO  35 
IFtCAOG  .GT.  f.5C.l  GO  TO  35 
FTCO=FT  (COO ) 

PTCPsFT I C  A 0° > 

FTCH=FT (C AOM1 

C.A1-CAO-?  .'0FL*FTC0/  (FTCR-PTO*) 

TP(A0S1FA11  .GT.  l.eri  GO  TO  35 
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CAlD=Cfll*0 
CA1P=A0S (C A1C  ) 

CAlG  =  A8«?(CAt*G) 

IF ( C AlO  .GT.  650.)  GO  TO  35 
IFCCA1G  .GT.  650.)  GO  TO  35 
FTC  1  =  FT  (C  A  1 ) 

IF  (AOS  (  (CAl-CAO)/CAO>  .LT.  .  J  0  1 )  GO  TO  10 
C  AO  =CA 1 
M=M  +  1 

IF(M  .GT  .10)  GO  TO  35 
GO  TO  5 

10  IF  (L  .FO.  ?)  GO  TO  20 
C  A=  C  A 1 

IF  (A9S(FTC1)  .GT.  1.)  GO  TO  35 
IF  (CA  .LT.  0.)  GO  TO  35 
FPCI)=FF 
CAP  < I )  =  CA 
1=1*1 
L  =  2 

GO  TO  30 
20  CA=CA1 

IF  (A°S(FTC1)  .GT.  1.)  GO  TO  35 

IF  (CA  .LT.  0 . )  GO  TO  35 

FM(J )=FF 

CAM( J)=CA 

J=J  +  1 

IF  (J  .EO.  T)  GO  TO  15 
IF  (j  .GT.  T)  GO  TO  31 
I  =  J-1 
K=K-1 
31  J=J-1 

15  K=KD 

IF  (K  .LF.  NF)  GC  TO  50 
PRINT  63 
11=1-1 
J1=J-1 
GO  TO  2  4 

35  PRINT  *  ,-ITFRATION  DOES  NOT  CONVERGE.  F=  ".EF,”  S  =  ~  , 
IF  <L  . EO .  2)  GO  TO  15 
L  =  ? 

GO  TO  ? 

2U  CONTINUE 

PRINT  63, (FP(T), CAP (I), 1=1, II, 10) 

PRINT  6-,<FM(J),CAM(J),J=l,Jl,10) 

PRINT  62 
DO  22  J=1 , J1 

IF  (J  .EQ.  1)  V0M(J)=5./PI*  (CAM(  2>-CAMl>  >/FO  EL 
IF  (J.FQ.  Jl)  VGH(J»=5./PI*fCAMCJl)-CAM(Jl-l))/F“£L 
TF  (J  .NF.  1  .AN0.  J  .NE.  Jl)  VGM(J)  = 

V  5. /PI*  (CAM(  J  +  D-CAH1  J-ll  >/FQEL*  .5 
22  CONTINUE 

PRINT  65, (FN( J) , VGH< J),J=1, Jl,10) 

PRINT  6n 


DO 

?1 

1=1, 

It 

IF 

(I 

.  FO. 

1  ) 

VGP  < 1 1  = 

5. /PI* (CAP ( 2)-CAP< 1) )/FOEL 

IF 

(  I 

.  FQ. 

ID 

VGP(I) 

=  5. /PI*  (CAP  (ID -CAP (1 1-1  D/FOEL 

IF 

(I 

.NE. 

1  . 

ANO.  I 

.NE.  ID  VG  F  ( I )  = 
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X5./PIM  C A  P  (T  ♦  1 )  -CAP(I-l)  >  /FD£L*.  5 
21  CONTINUE 

POINT  P  7  ,  (FP(I),VG»<I>, 1*1, 11,10) 

JQ= Jl/2 
CAM  0  =CA  M  (  J  0 1 
DO  23  J=1 , J1 
CAMN  =  CAMO-f?  M  f  J  ) 

CAMN=CIST*CAMN/PI 
C  AMN=  AM  00  (CAMN, 2  .  ) 

IF  (CAMN  .Gt.  1.)  CAMN=CAMN-2. 

IF  (CAMN  . lT.  -1.)  CAMN=CAMN*2. 

VNM  (  J)  =  lfl 0  CAMN 
23  CONTINUE 
PRINT  PC 

PRINT  66, (FM(J),VNM(J),J=1,J1,10) 

S-l  . 

00  25  I =1 ,Tl 
FF-FP (I ) 

CA=f AP( II 

0*2 . *°I *EP*1  .  F6 

PM-T  F/ ( 2.  «*1 7  t  0  •  > 

Ull=l.-0MH/( C** *2-OMH**2> 

U  2  ?  =  U 1 1 

U1  2  =  Of*/  (0M*»  ?-0MH**2  ) 
o=SQCT ( U1 1/U22) 

ALl=U22*fl+F  012 
AL2=U22*R-S‘U12 

PI  =  (T  <C  A )  ♦  i«  I  *  *  2*  (TANMCA*FLI-CA*EIM$ECH  (CA*EL  >  >  **2> 
P2=(Pl(rA)'oXp(Cfl»G)  -  C2(CA)*ExP(-CA*G))#j2/*», 

XMC0TF<CA*T1>  -CA*T1MCSCH(CA*T1>  )**2) 

P3=  .25*  <°1  <CA>**2)»  (EXP(2.*CA«G)-1.  >-.25MR2t  CA>  *‘2) 
XMFXP(-2.*CA*r>-l.)-Rl(CA)*P2<CA)*CA*G 
Pu=  AL1*  (T  f  C  C  ) '  *  2  )  *  (EXP  (2,*B''CA*D)-1.)-AL2*  (EXP(  -2 .*  R*CA*L)  -1.  ) 
X-u.*P»*2*U2?*r'A*0*T(CA) 

GEC  A  =  GE  (C  A  ) 

0P(I)=  0»U0  GECA  ++2* (Pl*P2*P3*P4)/o . /CA  **  2* -5 

PP  ( I ) =4  .  *  PC  (  T’ 

op (I )=ARS (FP ( I ) ) / ( <1 .-ETA) ♦  < 1.4ETA) *EN**2) *4. 

P.P(  I  )=F  ACT  •OP  » I ) 

26  CONTINUE 
0  =- 1 . 

no  2f  j  =  i,ji 

0P  =  PM(J) 

0=2 .*PI*EF*1. - 6 
CA  =  CAM<  J) 

OMr r  F/ (  2.  0  *1 7 r 0.  ) 

Ull= 1 OMH/(OM**  2-0MH**2 ) 

022-011 

U12  =  0M/ (OM**  7 -0MW**2 ) 
n=S0PT(Ul  1/022) 

At.  1=  U2?  *  p  ♦  F*  U  1  2 
AL2=U22*a-S*Ul? 

PI* IT (CA) ♦! . » *  *2* lTANH(CA*El> - C A  * E L * ( SECH  < C A*  EL )  )  *2) 

P2= (P 1 ( C A ) *rXF (C  A*G)-E2(CAI*cXP(-CA*G)>**  2/4. 

X* (CnTM  f  C A  *  T ! 1 -C A*  Tl* fCSCH (CA‘T1))**2> 

P?=  .  26*  (PI  fCA  ***  ?  )•  (EXP(2.*CA-G)  -1.  I-.2  5-  C  *?  2  (  C  A  »  *  2) 
v»  (rXO(-2, *  r  A  *  C ) - i . )  -  R 1 ( C  A  )  *  R  2 (CA)*CA*G 
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P4= ALi*  <T(CM**2»ME  XP(2.*B*CA*D1-i.  )  -A12*  <  EX  P< -2  .*  F*CA*D1 
X-4.  *8**2*l)22*CA*0*T<CAl 
GFC  A  =  GF  <C  A  ) 

PM(J1=  0*U0*  GECA  **2* ( Pi *P2*P3 *P4> /  . /CA** 2  .5 

PM(J1=4.*PM(J1 

RM< J»=ARS  <pm<  J>  >/ < <1 .-ETA) ♦( l.*ETA)  «EN**21  *4. 

PN< J>  =FACT*PN  <J> 

26  CONTINUE 
PPINT  60 

PRINT  71 ,  <FP( I ) ,RP(I )  ,1=1.1 1,10) 

C  PRINT  72, <FM(J1,RM(J1,J*1, Jl, 51 

IF  (II  .NE.  J 1 >  GO  TO  90 
OP  «<,  1  =  1,11 

84  PT(I )=RP( I)»PH(T) 

POINT  60 

PRINT  75, (FP(TI,FT<I1, 1=1, II. 101 
FP1=FP< 11 
FPl=FP< 111 

CALL  FTRAN(PT,PX,I1,FF1,FF11 
PRINT  60 
MrrNIM 
M=I1-1 

PRINT  73, (FP(I) ,PX<I1 ,1=1, M  ,201 
DO  54  1  =  1,  n 
XL=  2.*PI *FP ( T 1 *FLL 
XL=XL*1 .F6 
L  =  T 

IF  (I  . FQ  •  11  L  =  2 

SERP  (II  =  2  0.*  ALCG1Q ( <4.*RP< II /RG1/ 
v  (Cl.*  <RT  <T1  ♦RLI /RC-1**2 

X  ♦ C (PX<I1 ♦XL)/RG1**2) 1 

SERP <11 =SFPP(I1-76.4*DEIH*0IST/ 

X  (  2 .  *PI*  i.Ef  (FP(I^lt-FP(L-l)  >  /  (CAP  <1*11 -CAP  (L-l  M) 

SFRM(I1  =  2  0.*».LOG10<  <4.*RMm/RGl/ 

X  (  (1  .♦ (pt <11 ♦RL1/RG1 **2 

X  ♦  ( <PXIIH-XL!/RG>**21  1 

SFPM  (T1 =SEPN ( T1 -76.4*CELH*0I$T/ 

X (2." PI* 1 ,F6* <FM(I+H  -FN(L-l)) /(CAMdfll-CAM (L  - 1 1 1 • 

54  CONTINUE 
PRINT  60 

PRINT  77, (  F°  <  1 1  , SERP (I)  ,  I  =  1  ,  M  ,101 
PRINT  60 

PRINT  7*. <FP(I1,SERM(I1 ,I=1,N  ,10) 

TNT  N  =  35 on . 

OX=  ^0. 
yMTN  =  2500  . 

X  (*I  N=?4  0  0  • 

F  PR  GC=  .  01*  F°F  G 

XNIN  =  ATNT ( FRF  CO ) *  100  * 

X X  =  A  I N T ( . 01“FPI1-AINT (. 01*FLO) *1  . 

rx= i oo. 
ynin=o. 

ov=  2  or. . 

PY=  100  00  . 

no  c'.fe  j  =  i,ji 

IF  ( C AN ( J 1  .GT.  1  000  00.1  CAM( J)  =  100  000. 

IF  (VGM(J1  .GT.  lGOli.l  VGM ( J1  = 10  00 • 
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Pfe  CONTINUE 

DO  9 7  1=1,11 

IF(CAF(I)  ,GT  .  1  0  0  0  0  0  .)  C  AP  ( T  >  =  1  Ou  0  00  . 

IF  (VGP(I)  .GT.  1  000  .)  VGP ( I ) = 10  0 0 • 

P7  CONTINUE 

CALL  PLOT  (1.5  .Q.,-3) 

CALL  SYM«TOU  .  5  ,9.8,.  i.HEAO,  0,  20 
CALL  SYMBOL  < . 5 ,9.6,. 1  ,HEA01 , ) ,70 > 

CALL  SYM0OL  < . f ,c . 1 ,HEA02 ,0,70) 

CALL  SVMROL(.r ,9.?,. 1, HE ACT, 0,70) 

CALL  AX  IS < 0. , 0.  ,21HWAVE  NUMBER  (♦)  ( 1 /M) . 2 1 , 1 0 . , 9 1 . , Y N I N , D Y , 1 0  .  ) 

CALL  AXIS  (r  .  ,0.  tlGHFREQUENCY  (MHZ), -15, XX  , Q . , XM IN ,0 X , 1 U . > 

CALL  LTNE(Fm,CAM,J1,1,0,1,XHIN,DX,YMIN,DY,.08) 

CY= IOC. 

CALL  PLOT  (19.  ,0.  ,-T) 

CALL  A  X  I c (0. 

Y  ‘’5HGC>0UP  DELAY/CM  (♦)  (NSEC)  ,  25,  10  .  ,  .0  . ,  YMIN,  DY,  10  . ) 

CALL  AXIS  (n . , C . , 15HPPEQUENCY  (MHZ), -15, XX  , 0 . , X M  IN ,0 X , 1 0 .  > 

CALL  LINE  (PM, VGM,J1 ,1,0,1, X  MIN, DX, YMIN, DY,. 08  ) 

DY= 10000. 

CALL  PLOT  (19,  ,  0  •  »  -  T) 

CALL  AXIS (0. , r, ,21HWAVE  NUMBER  (-)  (1/M), 21,10. ,90.  ,YMIN,DY,10.) 

CALL  AXIS  (0.  ,0.  .15HFEECUFNCY  (MHZ), -15, XX  ,0  .  ,  XMIN  ,DX  ,  10  .  ) 

CALL  LINP(PP,CA°,I1, 1,0 ,1 ,XMIN,OX,YMIN,DY,.08) 
dy= ion. 

CALL  PLOT  (16.,  0.,-T> 

CALI  AXIS(C.,0., 

x  C5HGP0UP  OELAY/CM  (-)  (NSEC) ,25, 10.  » HO . , YMIN.DY,  10. ) 

CALL  AXIS  (  ;.  ,0.,15HFREQUENCY  (MHZ), -15, XX  , Q . , XMIN , D X , 1 3 . ) 

CALL  LINF (PP,VGP,I1,1,0,1,XMIN,DX,YMIN,DY,.C8 ) 

YMTN  =  -1  FOI  . 

OY=  36. 

CALL  PLOT  (IF  .  ,  3  . , -3) 

CALL  AXIS ( 0.  ,2?HNOPMAL  DISPERSION  (♦>  , 2 2, 1 0 .  ,  ' 0 . , YMI N , DY ,  1  0 . ) 
CALL  AX  T  S  (i.  ,0. , 15HFR£QUENCY  (MMZ),-15,XX  ,0  .  ,  XMIN »CX.  10  .  ) 

JJ=9 

58  CONTINUF 

DO  59  J=l,  J1 

IF  UJJ+J >  .GT.  Jl>  GO  TO  57 
FN( J )=FM( JJ+J ) 

VM( J)=VNM( jj* j) 

IF  ( J  . EO.  1)  GO  TO  59 

IF  !V“(J)  .  L  C .  VM(J-l))  GO  TO  59 

GO  t p  57 

59  CONTINUF 

5  7  J2= J- 1 

CALL  LINE  (rN, VM, J2, 1.0,1, XMIN, DX, YMIN, DY, .08) 

if  T  (JJ*1)  .LF,.  Jl)  GO  TO  55 
ymin=o. 

OY=  TO. 

00  27  1=1  .11 

IF  (  RF  f  I )  .GT.  300  .)  RP(I)  =  300. 

I F  ( R  P  ( I )  .GT.  200C.)  PP(I)=2000. 

?7  ~  ON  "■  T  NU  E 

CALL  Fi  OT  (19  .  ,  0  .,  -T) 

FALL  AYTS  C .  ,0  .,27MRAD.  °ES. .MINUS  WAVE  (OHMS) ,  27,  10 . ,90.  , 
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X  YMIN.OY.M.) 

CALL  AXTS  t".  .0.  ,15HFFECUFNCY  (MHZ), -15, XX  ,0  .  ,  X  M  TN  ,CX  ,  1  0  .  ) 

CALL  LINE (FP,PP,I1,1 , 0 , 1 , *M IN , OX , YM IN, 0 Y, .08) 

00  2  6  J  =  1 , J1 

IF  <FN<J>  .GT.  33  0. >  RM ( J )  =  33  3  • 

IF  ( P  M  { J  )  .GT.  2  C  0  0 .  )  PM  (  J ) =2000  • 

IF  <PT(J)  .GT.  300.)  RT(J)=33C. 

IF  (RT(J)  .GT.  2000.)  RT(J)=23QQ. 

26  CONTINUE 

CALL  PLOT  (17.  ,0.,-3) 

CALL  AXIS  (0. ,O.,?7HRA0.  RES..  PLUS  WAVE  (OHMS)  ,  2  7, 1 J  .  ,  90  .  , 

X  YMIN.OY.lb.) 

CALL  AXTS  (O. , 0. » 15H FREQUENCY  (MHZ), -15, XX  ,0 . , XMIN ,OX ,1 0 . ) 

CALL  LINE(FM,PM,J1,1,0,1,XMIN,DX  ,YM  IN  ,DY  ,  .  0  8) 

CALL  FLOT  (17.  ,G.,-3> 

CALL  AXIS  (  0.  ♦  C.  ,22HRAC.  RES.  TOTAL  ( CHMS)  ,22 , 1 0 .  ,  •  0 . , YHI N , D Y ,  1  0 . > 
CALL  AXIS  (0 . . 0 . . 15HFREQUENCY  (MHZ), -15, XX  , 0 . , X MIN , CX , 1 0 . ) 

CALL  L I  N£  (  FP  ,  c  T  ,  1 1 , 1  ,  Q  ,  1  ,  XM  IN  ,  DX  ,  YM  IN  ,D  Y  «  •  0  8 ) 

I1=M 
J1  =  M 

YMT  N  =  -l 0  O  0  0 . 

YMI N  =  -2  50  . 

OY=20C0  . 

DY  =  5  0  . 

DO  P2  1=1  ,  Tl 

IF  (PX(T)  .IT.  -13000.)  PX(I) s-10000. 

IF  (PX(I)  .LT.  -250.)  PX  (I ) =-250  • 

IF  (PX(I)  .GT.  13)00.)  PX(I)=10000. 

IF  (PX(I)  .GT.  250.)  PX  tl )  =  2  5 0  » 

9?  CONTINUE 

CALL  FLOT(17.,0.,-3) 

call  axis ( 0. , C . ,22HPAC.  RFAC  TOTAL  (OHMS)  ,  2  2,  1  0  .  , no  . ,  YMI N  ,  OY  ,  1  0 .  > 
CALL  AXIS  (0.  ,0.,15HFR£QUENCY  (MHZ), -15, XX  ,0  . , X  MIN ,DX , 1 0 . ) 

CALL  LINE(FP,FX,I1,1,C,1,XMIN,DX  ,YMIN,DY,  .0  6) 

on  03  1=1,11 

IF  (SE  PP  (  I )  .LT.  -#»0.  )  SER°  ( I  >  *-8  0  » 

IE  (SERM(I)  .LT.  -60.  )  SERM(I)=-e0. 

93  CONTINUE 
YMTN=-8T. 

O  Y=  10  . 

CALL  PLOT  (ie.  ,0.,-?) 

CALL  SYMOOL(.r,4.8,« 1, HE AC, 3,20) 

CALI  SYMPOl ( . e ,9.0,. l.HEAOl ,3, 70  ) 

CALL  SYMPOU  .5,9.4,.  l,HEAO2,3,70> 

CALL  GVMOOL( . 5 ,Q. 2 ,. 1  , HE A 03  ,0 ,70  ) 

CALL  AXTS (0. , f. .26H-INS.  LOSS, MINUS  WAVE  (OB) ,2b,  e.,90.,YMIN, 

X  OY  ,10  .) 

CALL  AXIS  (’ .  ,3.  .15HFREQUENCY  (MHZ), -15, XX  , 0 . , XM IN ,0 X, 10 . ) 

CALL  LINE (FP.SERP.Il , 1,0 ,1,XMTN,CX, YMIN.OY, .08 ) 

CALL  PLOT  (IF.  ,0.  , -3) 

CALI  SYMBOL  (  .  r  ,  9. 8,.  1  ,HEAn  ,  £1,23) 

CALL  <rYM°OL(.5  ,9,f  ,.  1,  HE  ADI, 3, 70  ) 

CALL  SYMBOL  (  .5,9.4,  .  1 ,  HEAD  2 , 0 , 70  ) 

CALL  SYMBnL(.5,9.2,.l,HEA03,0,70) 

CALL  AXTS  (0. . ? .  ,26H- INS.  LOSS,  PLUS  WAVE  (OB)  ,2b,  6.,S0.,YMIN, 

X  OY.IO.) 

CALL  AXIS  (n.  ,0.,15H FREQUENCY  (MHZ), -15, XX  ,0 . , X M IN , C X , 1 0 . ) 
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CALL  LINS  (FM.'TRM.Il  ,  1 , 0 , 1 , XMIN , OX , YMIN , OY ,  .0  8) 

TAIL  ENOPLT 

stop 

<50  P°I NT  76 
STOP 

60  FCPMAT(lHi) 

61  fn?M«T(5X,"  Hr  **,  £15 •  7/5X  ,**  Tl  =  **,  £  15  .7/*;  X  t  -  C=  *',E15 .7/ 

X5X,"  G=  ** «  El5 . 7/ 5X  »**  L  =  ",E15.7/ 

X  5  X ,  **  N=  ",£1F.7/5X,"E7A=~,E15 . 7////9X  ,“NO.  OF  F'S  APE  M,I5) 

62  fopmAT(//1C<F,:15.5/)  ) 

63  FORMAT(///5"X,**S  =  l"// <10X,“F=  ** ,  E  15  .  ?  ,  1  OX  ,  *‘K  (- >  -  «,£i 5.7/)) 

64  FORMAT(///5''X,**S  =  -1**//(10X,**p  =  **  ,  E15 . 7 , 1 .  X  ,  "K  ( ♦ )  =  ",E 15.7/)) 

66  FO»MATf /// »10X,**F=  ",£15.  7,13X,**F  (~>  =  **,  El  5.  7/ ) ) 

67  POPMAT(///  <1CX,"F=  ",El5.7.10X,**P  (  +  3  =~,El5.7/)> 

66  FOPMAT(///?y,~li=  **,E15.7/5X,**A  =  ", E 15. 7/5X ,“P=  ",El5.?/ 
X5X.**TC-  -,F15  .7/5X,**N=  **,  E 1 5 .  ?/ 5  X  ,**ET  As  ** ,  E 15  . 7  ) 

71  FCFMAT(///(10V,*‘F-  ** ,  E 15 . 7 , 10 X  ,"R  AO .  RES.  <->  =  ",E15.7/>> 

72  F0=>MAT(///(1CX,"F=  “  ,  E 1 5  .  ?  ,  1 0  X  R  AD  .  FES.  (♦)  =  **,El5.7/>> 

73  FORMAT!  //flOX,*’F=  ",Ei5»7,10X»" RAO.  REAC  *  TO  T  A  l  =  "  ,  E  15 .  7/ )  ) 

75  F  ORM  AT  (//  /  <10X,*‘E=  ",E15,7,1JX,‘*RA0.  RES.  TOTAL  =”,El5.7/>> 

76  ffRMAT  (**l  A  X  ROOT  EXISTS  ^CR  CNE  NAVE  ONLY**) 

77  FOPMATt// (ir.y,*«f=  ", El5. 7, 1QX ,"I*S.  LCSS  (-)  =  M,E15.7/)) 

78  fOPMAT<//aaX,“F=  **,E15. 7,lCX,"ItsS.  LCSS  ( * )  =  “,E15.7/>> 

«1  F  0s  M  A  T  I  /  /  (i^X  , "  L 1 1" » I  ** ,  ”  3  =  " »  El5 • 7/ ) ) 

8?  FORMAT* //T1CX,"  At*M4,“>  =  ",E15.7/>> 

63  FOOMATT// UOX,"  PI",  14,**>=  ",El5.7/>) 

85  FCP  M  A  T  I  /  I  IPX  F=  **,E15.7,10X  ."GROUP  CELAY  (♦  )  =  “,E15.7/>> 

86  FORMAT!"  r=  ",E15.7,"  NORM  DISPERSION  <♦)  =  **,F15.7» 

87  FOP  MAT  I  **  F  -z  "  ,  F15  •  7  GROUP  OELAY  <->  =  **  *  -S 1 5 . 7  > 

102  FOPMATC2A10) 

100  FQPMAT(7A10I 

101  FORMAT ( 7  A  1 0) 

END 
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FUNCTION  SECH(CA) 

COMMCN  FL,AL1,AL2,3,C»T1»G,S*ETA»EN»P,AY,A 
SrCH=C. 

IF  (CA  .L  E .  7 1, 0  . >  SECM1  ./COSH  (CA) 

RFTUPN 
r  MP 


FUNCTION  r SCH ( CA) 

COMMON  EL  .All ,AL2»8,0»Tl ,G,S»£TA,EN*P»AY»A 

csch=c. 

IF  (CA  .LF.  7 A 0 . >  CSCH  =  1./SINH (CA) 

KFTUFN 

FNO 


FUNFTTON  rnTH(CA) 

COMMCN  FL , AL1.AL2 ,3,0.Tl,GfS, ETA  ,EN  ,P«AY, A 
COTH  =i./TANH(CA) 

RCTUPN 

END 


FUNCTION  T(CA» 

COMMON  Fl.ALI ,AL2 ,T,0,T1,G,S,ETA,EN,P,AY. A 
T  =( AL?+T ANH (CA*EL) )/ (AL1-TANH(CA*EL)) 
PFTUPN 
END 


FUNCTION  Rl(fA) 

TOMMON  FL , Al 1 , AL2 ,3,0 ,T1 ,G,S, ETA  ,EN  ,P  ,AY , A 
Rl  =( l.-AL?)*EXP(-B*CA*0» ♦ (1 . ♦ A LI >* T (C A ) * E X P ( b* C A*0> 

PFTURN 
FNO 


FUNCTION  P?rCA) 

COMMON  FL, ALi, AL2,3,0,T1,G, S, ETA,EN,P,AY,  a 

R?  =  (1  .♦AL?)»EXP(-n*CA*(m  (l.-ALl»*T(CA)*EXP(e  C A*CI 

RFTUPN 

FNO 


FUNCTION  FT(Cfl) 

COMMON  EL  »ALl,AL2,O,0»Tl,G,S»ETA  ,EN,P,AY,A 
FT  =  ,F«  ( (COTH  (CA*Tl )-l . >  *R2 (C A)*EXF  (-CA*G> *E  X  F ( -P*  CA*  D) 
X-  (COTH<  CA-»Tl )  »1 .) *R1 (CA)  *EXP(CA*G)*EXP(-P*CA*0)  ) 

RCTUPN 


FUNCTION  FT1 ?rA> 

COMMON  FL»ALl*AL2»8,0.TitG»S»ETA*EN»P»AY.A 
f XT  -p*p( _f*cA*0>  * (R1 CCA) *EXP  «CA*G)-R2<CA»  *EXP  <-CA*G> ) 

X  *  S*  T  1  *  (OSCMfCA^T  1 ) )  *  *  2 


x.s.r,rxP(-o*  CA*D>  *  <<CCTH(CA*Ti)  +  i.  >*R1  (CA  )*EXP  (CA  G) 
tflCOTMU‘UI-1.  >*R2  CCA)  *EXP(-CA*G)  ) 

FT  T  T  = 

¥  4.3.  *o»s*0,,'rXP  (-2.  *3*CA*C)*  (  (COTHCCA*Ti)  ♦!.)*(  l.-AL  2) 
X*  £YP  tCA‘G)  -  (COTMCA*T  1)  -  i. )  *( !.♦ AL2)*EXP(-CA*G>  > 


X>S*T  CM  AL1  +  ALI1* (SECMCCA*FL)**2>/ (AL1-TANH(CA*EL>  >**2 
FTT  F  = 

X  CCrnTMCCA*Tl »-l. »*  Ci.-AL1)*EXP (-CA*G>- (COT H C C A* T 1 ) *1 . ) 
X*  (l  .♦ALD  *  FYP  (CA*G)  ) 

FTT4=FTT^*»"rT  F 

FTI-.5» (FTT+FTT24FTT3+FTTA1 

PFTUPN 

ENO 


FUNCTION  CINC(CA) 

COMMON  -L,ALl,AL2.8,O.Tl»G,S»ETA,EN,P,AY,A.ELl(4C>.PE(4i)>,AA(40) 
PI*  5.1%  1  “502054 
$  INC  =  <F  IN  ( PT  *  0  A  )  )  /(P  I*CA) 

RETURN 

FNO 


FUNCTION  GAY(CA) 

COMPLEX  C.CF 

FL*flL1’aL-*g*O’T1,G*S’ETft»EN»p»AY*A.ELl(40>,PE(40),AAU0> 

T  *  L  “  □  0  r 

PI=T. 141502654 

n=fn 

r=CMPLx  (o . ,  "i . » 

00  1  1=1. N 

CS=CMPL  X(OOc'(CA*I»PE(I)),-SIN(CA*I*PE(I))> 

IF  ( LMOOr  .FO.  2)  GO  TO  2 

C=r«-SINC(  .  5'  AA  (I)*CA/FI)*FT  A**I*SQRT  (ELK  I)  )*CS 

GO  TO  5 

2  CONTINUE 

C=045INC(2.AA(T)/(PE(I)*(3.-ETA)))* 

XSINC <<CA* PF< I* «.5/PI)  -.25  * (3  ,4ETA))*FTA** I*SQc T  (EL1(I) )*C  S 

3  continue  “ 

l  continue 
G AVr c APS  f  r > 

RETURN 

ENO 


FUNCTION  GF(CA) 

COMMON  FL , Alt , AL2 .8,0,T1 «G. S, ETA  ,EN  ,P,AY, A.EL 1 (40) ,PE(40) »AA (40) 
GF  =APG(GAY(CA)  *EXP(-3*CA4D)  /FTl  (CA)) 

RFTIJPN 

fnt 
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o<MT 


SU3R0UTINF  H TRAN ( R , X , N , F 9EG ,F END ) 

DIMENSION  R(3),X(3) 

PI  =  3.i4i51=2f'c?5S 
Fn£l  =  (F"ND-FFPG)/  (N-D 
F  =  F3FG«- .5*  FPFt 
INC=MOO (N , 2) 

NI=  N  +  INC-  I 
NM 1 =N- 1 
NIM?=NI -? 

00  33  I=1,NM1 

xm=o. 

IF  (I  ,FQ.  1)  ®  X  = ( 3  •  *  R (1) ♦  6  •  *  R (  2  )  -  R ( 3 ) ) / ft • 

IF  (I  .FQ.  NM1)  RX=  { -R  (N-2) +6.*R  CNHD  *3.*F  (N)>/8 
IF  (I  .  EQ .  1  .Oft.  I  .EQ.  NMD  GO  TO  20 
RX=  <-R(I-D*9  .*R(I)+9.*R  (I«-D-R(I*2)  »/16. 

CONTINUF 
FI  =  F  °FG 

00  26  I P=  1 , NT  N  ? , 2 

X(D=X(  ri  ♦4,*(R(IP*l)-RX>/( I  F I+FDEl  )**?-F**2  ) 

X  ♦2.*{R(IP  )-RX)/<  FI  **2-F**2) 

F I-F 1 *2 •* FDrL 
CONTINUF 
F£N=  FENO 

IF ( I NC  .EO.  C  )  FEN  =  FENC-FOEL 
X (I )=v< I) ♦ (p  < M >-RX)/ (FEN** 2-F**  2) 

X  -<®(1  )  -RX) / (F6EG**2-F**2) 

X(T)=FDFl/3.*>(I> 

TF(TNC  . FO  .  1 1  GO  TO  30 

X  (I ) =X( I) ♦,D*FOFL* ( (R (Nil -RX) / (F EN** 2-F**  2 ) 

X  ♦  <c  (N)-RX)/ <FEN0**2-F**2)  ) 

X( 1)=?.  /P  T*F*  X  (  I )  +  RX  /PI*  A  LOG 
X  <  < 1 ,-F/FENO  > /  (IDF/ FENO)  *<F«-F8EG>/  (F-FQFG)  ) 

F -F ♦FDEI 
rcNTiNur 

NM  ?  =  N -? 

Xl-(15.*X(l)-10.*X(2)+3.*X(3))/6. 

X2=(3.‘X(D+F.*X(2)-X(3))/«. 

00  31  1=3, NM  ? 

X  J-  <-y  (1-2)  4-9  .*X  (1-1  )  +  9.*X  (I) -X  (1*1)1/16. 

X (I-?)=  XI 
XDX2 
X2=XT 
CONTINUF 

X (N )  =  (1 5. * X ( N M 1 )-10. *X  <NM2)+3.*X (N-3) >/8. 

X  (N  “  1  )=  (3.*V(NM1D6.*X(NM2)  -X(N-3))/8. 
x (N-2)=X2 
X  (N  -  3  )  =  X  1 
orri/PN 


B.  Microstrip  Model 

A  second  computer  program  incorporating  the  microstrip  model  has 
also  been  completed  for  the  CDC  6600  at  Hanscom  AFB,  Ma.  The  physical 
quantities  graphically  displayed  by  this  program  are  the  wave  number, 
group  delay  and  insertion  loss  for  both  solutions,  normalized  dispersion 
for  the  +  wave,  input  resistance,  corresponding  reactance  and  the 
magnitude  of  the  impedance.  Print  out  is  provided  by  the  program  as 
for  the  basic  theory  model. 

Note  that  apodization  and  normal  modes  are  not  permitted  here  and 
that  additional  input  constants  are  required.  Thus  the  use  of  the  input 
cards  are  modified  from  the  basic  theory  program  by  the  following: 

Cards  2-4  -  The  increment  values  are  always  o  and  the  option  values 
should  be  o. 

Card  5A  -  Zcl,  Fc,  o,  acK 

This  is  a  new  card  to  be  inserted  between  card  5  and  card  6.  The 
indicated  constants,  required  by  the  microstrip  model,  are  to  be 
inserted  here ,  separated  by  commas . 

Card  6  -  columns  1-70  may  be  used.  The  first  ten  columns  should 
contain  IND  COND_. 

The  listing  of  the  entire  program,  except  for  control  cards  which  are 
the  same  as  for  the  basic  theory  program,  now  follows. 
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v :  -  y  ■ 
y  y ’ 

7  1  '  ''NT  T  MJ 

Y  ( N  )  -  ♦  1  -  .  lX(NM1  >-n.*X(NM?»  *  *  .  *X  (N-  3  )  )/C. 
vn-:i:(7.*y(M;)*6,*)((NM?)-X(N-3))/T. 

y  ('■  -  o  )  -  v  ? 

Y  (''-’hVI 
F  “r  I  I  r  N 

f  '  “ 


F|  rN'‘  T  TFI‘!  I  C'  ) 

r  PMwrv  rL  *SI  l.Al2,X,D*Tl,Gt5fETAtENtp»AY,AtEl  l('*0)»FE(40)»Afl(«r)) 
g  =aff  (r,flv(rfl,  *  fxpc-x+ca-oj/ft:  (Cfi>> 

•3  r  T  I  ■  D  K 
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C.  Generalized  Dispersion  Relation 

Another  program  has  been  implemented  on  the  CDC  6600  at  Hanscom 
AFB,  Ma  which  solves  the  dispersion  relation  problem  for  surface  waves 
including  any  arbitrary  orientation  of  the  biasing  field.  Computer 
print-outs  and  plots  of  wave  number,  wave  length  and  group  delay,  as 
functions  of  frequency,  for  both  +  and  -  solutions,  are  provided  by 
the  program. 

The  input  cards  to  this  program  are: 

Card  1  -  Hq,  t^,  d,  0,  £,  <p 

These  six  quantities,  separated  by  commas,  are  supplied  here.  The 
lengths  are  in  meters  and  the  angles  are  in  degrees. 

Card  2  -  first  f ,  Af ,  number  of  frequency  values 
Here  the  user  is  to  provide  the  first  frequency  value,  the  frequency 
increment  and  the  number  of  frequency  values  to  be  considered,  all 
separated  by  commas.  The  maximum  number  of  frequency  values  per¬ 
mitted  is  currently  500.  Although  the  program  could  compute  the 
frequency  range  itself  as  in  the  other  programs  it  was  left  for  input 
to  provide  flexibility. 

Card  5  -  heading 
Card  4  -  heading 

Two  heading  cards  for  the  graphs  are  here  required.  Columns  1-70  may 
be  used. 

The  listing  of  the  entire  program,  excluding  the  standard  control 
cards,  now  follows. 
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oht-o  orG  THrTf=c>U  DFG  - -  SURFACE  WAV'S 

unqr.cm  >/0L*<  TNRUT,  OUTPUT) 
p  T  m  r  K  c  J  ON  P-CGIt.  (3) 

^ion  F  ( 5  1 )  ,  FP  (  5-  1 )  ,FH  (SHI)  ,CAP(  5GD  ,  C  AM  (  5  0  1) 
OTM'-'NST  ON  H~  A  r  1  (  7  )  ,HE  AD?  (  7  ) 

X  ,VGP  (5H  1 )  ,VGH  (501) , ELAP (501) »£L AM (5  01) 

COMI-CK*  r,  -  L,Ti.  ,S,r,ALl,AL2 
r!  nf>  »,M,?l.r,THET,FL,PHI 
“•tar.  «,rn.r.,rr£L,NF 
fftC  l  o  n , h f a  r 1 

r  r  A"  |Q-  ,MFdT  3 

pc-[  MT  f  r  . 

P  -  TNT  f 1, H ,-l .P.ThET , EL , PHI 
P  r  I  N’T  t:  H 

p  :  t  r.  PFOGTC-'  FHEFTH4RES*AH2567tlGhJ«WEINPFFG/ 

CALL  FI.TTP  XtPROriO,  100  .  ,  12  .  ,  1  .) 

-  o  u  n  T  =  1  ,  Nr 
-  F  (  T  )  =r  PrT.t  (T-l  ) -FDEL 

5  . 

Tw.-;.th!-T.'  pT/1fQ. 

dmj . pf r ip i /  i  f  c . 

c  m  ’  7  r , 
r  AM  =  2.  - 

It-  r-AP*  '  ?  1  H  "M 

Y  T  « 

*»■  —  A 

J-1 

0  FFrP(K) 

THr (GAM‘H)*  2-FF**2 
rr 

i:vy  r  7 )  /  TQ-  (  .  I  N  (  T  H  f T  )  |**2»1. 

uyy  r  tl/  Tp  (  <rTN  (THET  )  *SIN  (PHI)  )  **2  +  (COS  CHET)  )  **2>+  1. 
uxv  V  vr-  ->.  •  T,t/Tn«cnS(  TFET)*COS  (PHI  )♦  STMTHF.T  ) 
ijr  1  v  VT;  -  ?/  T')'  r  IN  (THt  T  )  *  S  IN  (  PHI ) 

T'Hr||YYvX  **'■“>..  *'JXX*UYY 

T,'i^r-TCM 

l  -? 

IF  (TFM  .IF.  r.)  PPINT  NEGATIVE  SQUARE  R  GO  T  *' 

IF  <tm  .i  f,  ,  ,  )  r,n  TO  3  5 

i'rS'’CT(TrM)  .-/UYY 

l-l 

:t  IF  (L  .  rO.  ' )  r,n  to  2 
1  r  -  1  . 

GO  TO  3 

T  r  PM TTM|  ' 

Mt_«-r*||YV*E  I :  v  y  Y  Y  I 
'  1  -  p  -  1 1  y  v  -  C  L  T  X  Y  Y  I 

T  f  ( J  . G -  .  l )  GOTO  53 
•  n-  ( a  I  ■’♦1  .  )/  (  At  2-1.) 

PJO  -  ft  LOP-  (  t  Al  2-1  . )  /  (TC"  (AL1.  +  1 .  )  >  )  /2.  /C/O 
t i  (rrn  .  i  T.  .  )  CAO=  ioo  . 

3  rrNTmr 

IF  ( T  ,'n.  1  )  GO  TO  51 
I-  U  ,rG  II  GO  1C  51 
*FQ  .El.  1'  -■  AOrCAP  (1-1 ) 

:r(i  .fp.  2)  f ao=gah  ( j-i ) 
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-  ” r«.  = ,  i?  *r.  on 

(  Ano=r  a  o  +  nn 
r  aoM=can-oct 
r a^r' -ta n*  o-*c 
r  t  (r  *rr  > 

IP  (ram  ,r,T.  .-5-. .)  GC  7C  35 

PTCO=FT  (Gael 
r  Tp  o  =  rr  (C.  n  OP  ) 

PT'-^rT  (C  ACM) 

rAi=GfO-?  .  '  'CI'FTCO/  ( PTC P-PT Cy ) 

IPflopiGM)  .GT.  1.F7)  GO  TO  35 
r  a  I  r  -Cfl  i*r*r 
r  r  i  p  -  a  o  r  *  c  a « n  » 

TC( caio  . G T .  -GC.)  GO  TO  3 5 
rrr  i  =ct  (C.  f.l ) 

c  pr'T\,r  *  .rap.  PTCO.cai  ,  FTC1 

Tc  (AnG ( ( 0 £1 -PflO) /CA 0  )  .LT.  .001)  GO  TO  10 
raor  rai 

W  +  « 

T p  (  M  , G  T ,  1 o >  rn  TO  7K 
on  t  p  r. 

10  IP  (l  . eO .  °)  GC  TO  ?0 

CfliCAi 

T  p  CarG(PTr*)  .  G  T  ,  1.)  GO  TO  IP 
t r  <ra  . lt .  n  )  goto  35 
pp  ri ) =cr 
r.!\*  (r  jrca 
T  -  T  + 1 

:  p'Int  »  .-cA-  ",ra,"p=  ••  t ec 

i  -? 

GO  Tp 

?u  g  a = r  f i 

I c  1A=G(PTr-)  .GT,  1.)  GO  TO  35 

IP  ( C  a  .IT,  -’.)  Go  TO  3  5 

P>(  J )=r r 

r£M  ( j)=CA 

J-.J  +  1 

C  p " T k T  *  ,,,Cfl=  *’,ca,"p=  •* , Ec 
IP  fj  .PO.  n  go  To  15 
IP  ( J  . GT .  I)  GO  TO  31 
T  =  )-l 
K  =<  - 1 
31  J-J-1 

15  K  r  <  4  J 

JP  ( K  ,  Gc.  MF '  GO  TO  50 

o?TMf  £ n 

T 1- T-1 

j:  =  j-  i 

GO  TO 

7,3  d  *j  I  m  t  ♦  ,”TT"  *3  AT  ICN  O  C£S  NOT  CCNVERGF  .  F=  "tEF,"  S=  **,S 
TP  C L  .50.  ?)  Gr'  TO  15 
i-  n 

Gr  JO  ? 

?U  r  'Mf  I  N'U  r 

prjT'T  p  ■»,  f  P”  (!  L)  .C1P  »LL»  »LL  =  1  .11  ,5) 

^Tr  (  U  ,  (PM(U_)  fcaM  (LU  »LL  =  1,J1,5) 
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sasK 


* 


pp 

1  1 

T-1 

'  t 

1  F 

f  1  )  - 

r  n 

L 

J-l 

A  y  ( 

J)  -2 

p - 

I  MT 

6  f  . 

p: 

IN  T 

p  p  . 

-P 

r  ■* 

T 

JF 

<  T 

.  fp 

I  F 

<  T 

.  -n 

IF 

(  7 

.  ME 

• 

/  PT 

( G  A 

)  V  P-p  (  T  )  =  0 .  /  P  T*  (CAP(  2 )  -  C  A  F ( 1 >  )/FOEL 


IF  (T  .ME.  :  .  ANP.  I  ,NE.  Tl)  VGF<I>= 

.  /  PT  ( I';tp  (T+  ^  )  -CA'M  1-1 )  >  /FOEL*.  5 
'  PNT  T  NU  - 

P-  tn;t  >  7  %  (FP  M  )  ,  VGP<  I  )  ,  1=1,  II  ,  1(1  ) 

~ n  n  :■  j  =  i ,  j 4 

IF  (j  ,-r,  i \  VCM(J)=E./FI* (CAM( Z)-CAK(l) )/FDEL 
T  (J  .FT.  m  VGM(J  )  =5./PI*  (GAM  (J1  )-CAM(  Jl-1  )  ) /FPEL 

P  (J  .  MF  .  .<  .AMD.  J  .NE.  Jl)  VGMCJ)= 

./PI-  (f  AM  J  FI  )  -  CAM  {  J-l  )  )/Frj£L*.5 

~%r.  r  Mu  " 

T  r.T  «  r  ,  (F*’  (  J)  ,  V  r,M  (  J  >  ,  J  =  1 ,  JL  ,  1  U  ) 

N't  ■">  'i  . 

y- 1  nr. 

y  MT  M  r  o  1  m  . 

"  Y  -  '  > 

YMT'.|Tr, 

'  v  =  ■>  r : :  . . 

r  n  i  7  I  -  1 .  T  i 

re  »rar»(->  . g’.  y)na:.)  cap<i)=50ocgc. 

T  r  ( *■  l  5  p  (  T  )  .ft.  mno.  >  ELAF  (I  >  =  1G0C. 
f:  d  ( i )  -  '  r («r  f  ( I )  ) 

:r  (vr.T<T)  ,  r,  ’ ,  i  j  i  c. )  vg°(I)sico:. 

:  0  ;»  -j-l  .  J  - 

T  f  (  f  A  M  ( J  )  .  F  T  .  5  j'Lflii.)  C  A  M  {  O  )  =  5  0  0  o  u  0 . 

IP  riMI  n  ,GT.  1  ;  i  0  •  )  E  L  AM  J  )=1ujQ. 

>/rv(  (  =  (>/CM  (J)  ) 

TP  f  V  E  M ( J '  . G  T .  1  ICC.)  VGMfJ)  =  lU0o. 

v-n. 

1  ;  . 


r  "l  L  rVMP  OL  (  .  "  , c  1 , HF  AU? ,  1 , 70  ) 

t  U  ?YTSfr.»  .  ,  2  E  H  W  A  V  F  NUMFEK  KM)  <  1/MfT  EF  )  ,  2  6  ,  v  ,  90  .  ,  Y  MI  N  ,  CY  ,  YY) 

f  .'LL  iTTMi.  ,  .  » 1 5>H  F  P  OU  E  NC  Y  (MHZ), -15, x  ,  C .  »  XM  IN, OX , YY ) 

C  r  L  l  L  T  Nr  (Fy.CM, 

Jl,l,n,l,XMlN,OX,YNIN,OV,  .05) 

PAU  ri  n-  ( ;  t . .  i  . 

'•■'ll  s  v  >•  =  o  L  < .  ,  <  .  t- , .  i ,  nr  A  p  1  , ) ,  ?  g  ) 

r  r  L  L  GY  VP  01  (  .  •'  ,  -  .  4  , .  1  ,H£  AH?  ,  I  ,  7u  ) 

-all  AXIMC.,  , , 2  EHW  A  V  F  NUM  nFR  K  (  - »  (  1/ Mf  TE  P)  ,  2  6  ,  *,  50  .  ,  YM  IN  ,  0  Y  ,  YY  l 


vy  =  i 

p  . 

r  u  i 

FL  PT  M  . 

pal  1 

<-  Y  M n  |  ( 

r  "LL 

r Y MP  OL  f 

f  U  t. 

:  y  t  s  r r . 

r  ■'  L  L 

AY  T  P  (  .1 . 

f  •'  L  t 

LTMF  (P 

PAL  l 

rt  n*  ( ; 

'’LL 

5  y  p  l  f 

PALI 

TY  YPOI  < 

‘ALL 

A  Y  J  p  (P. 

p  1 L  I 

r  v  t  r,  (  p  . 

'  A  L  L 

L  IMF  (Ff1 

1 1, 1  ,  1,1 , XMTN ,0X  ,YMIN,OY. . jf) 


FPL!  rl  nr 

(ALL  FYM-OL  <  . -  ,5  .o, .  1  ,  HFA01  ,  )  , 70  ) 


******>  V  ■****& 


’  m+j  , 


CAM  'YM»  f'L  (  .  - ,  .  .4,  .  1  ,HF  ftn?  ,  ) ,  70  ) 

CALL  A  X  TS ( 0 .  /  iTi  HWA  V  ?  L  L  Is  C  T  t-  (♦  >  (  P  ]  C  C  C  N '  )  ,  24  ,  X  X  ,  0  .  ,  Y  M I N ,  0  Y  ,  V  Y  ) 

ff-ll  /X  IS  <  .  .  ,  n  1?HP“ECL-Nf.Y  (HHZ)  ,-15  ,XX  ,  0  .  ,  X  f  I  N  ,  D  X  ,  Y  V  > 

'"ALL  LIN"  (F  v.tfL  AM,  ji  ,  i,  o  ,1  ,X*IN  ,CX  ,YMIN,  IV  Q8  ) 

r all  pl  ot  <<  x .  ,  j.  ) 

nLl  cYM*ni  (  .  ’  ,E  .h,  .  1 ,  HF  A?!  ,  1 ,  70  > 

CALL  "VMOOL  t  .  ,-.u,  .  1  ,  HEAD?  ,  1 , 70  ) 

CALL  fix  IS  (  P.  ,  "  ,  ?t*  pWf  VFLtNGTH  (-)  (  M  1  C  P  ON*:  )  ,  Zu,  ,  X  X  ,  CO  .  ,  YMI  N  ,  0  X  ,  Y  Y  ) 

(.ML  fi  X  T  c  f .  n  .  ,  1FH  FrLCUENCX  (MF7),-15,XX  ,  0  .  ,  X  P  TN  ,  D  X  ,  Y  Y  ) 

r  AL  l  LIN"  (fP  .  CL  AP,I  1  , 1  ,  Q  ,1  ,  XMN  ,C>  ,  YFIN.  L  Y  , .  08  ) 
r 2 l i  pl  ("T  t  ^  . ,  o . ,  -X) 

("ALL  *"  Y  M  C  0 1  (  .  ~  .  p  .  F  ,  .  1  ,  HF  A  D 1  ,  )  ,  7  L  ) 

"!Ll  Cvmf>01  f  .  r  ,L  .  .  1 ,  HEAD?  ,  1 , 70  > 

c a l  i  f.v  if  f , 

X  "3HGFOUP  tJE  L  a v /C  N  (♦>  (NSEC)  *  25*  XX  ,  CQ  .  ,  YMIN  ,  1 Y  ,  Y  V  ) 

r  ALL  ?vrr  f  ■  ,  r  .  ,  1  5F  F  P E  C  Li c  N Y  <MFZ),-15,XX  ,  0  .  ,  X fTN , 0 X , Y Y ) 

"fill  L  I  Nc  f  t-M  ,vr,M  ,.J1  ,  1  ,  0  ,  1  ,y  f-T^  tux  ,  YMN,  CY  ,  .Cd  ) 

fit’  PLOTM 

^  L  1  CY  v‘"  rlL  (.  1  1  , H F  A 1  ,  1,70) 

r  r,  L  V  <"YM?nLf.T,  .u, .  1  ,  ME  An'5 ,  >  ,  7C  ) 

"'LL  A  X 1  -  (  .  ,  .  , 

x  ,  5HG  pOU  c  CE  L  A  Y  /  "  M  (-)  (NSEC)  ,  2  5  ,  X  X  ,  0  .  ,  Y  MI  N  ,  U  Y  ,  v  y  ) 

b.'U  fix  T  f  .  ,  0  .  ,  15HFCECUENC  Y  (MHZ), -15, XX  ,  0  . ,  X  Pi  N  ,  0  X  ,  Y  Y  ) 

r"L  l  i.r  ‘IT  (f-n,  vr,p,  n,  1 ,  u,  1 ,  P)  ,D  X  ,  YM  SOY  , .  id  ) 

I  ALL  "  N  n  p  L T 
r  -  qo 

F n  rr?MAT( 1K1) 

FI  Fp:  MflTf  ->Y  , "  p:  FlL  .  7/5X  Tl=  " ,  E 1 5 . 7/ "  X ,  "  C=  M,El5.7/ 

v  twcta  ,  Fir  .  y/3X,"  L  =  **,ri3  .  7/5X  ,"PH  =  ", ElE.fi) 

6.?  rr:wn(//,r(:rlr.5/)) 

F7  r  np  YfiT ( ///rnx , "S  =  1 " / / (10X,"F=  ",E15.7,inx,”K  (-)-  "tfcl5.7/>> 

6*  rA^HTI/// ^  y  ,"r  =-l"// (lux,"e=  "  ,  El  5 . 7 , 1  ;  X  ,  (♦)  =  M,El5.7/)) 

h  ff  F''-  MfiTt ///c  y."S=l"//(lOX,"F=  ",E15.7,1CX  ,"LAM(-)-  M,E15.7/>) 

Ff,  Fp-Mr(///r  X,**f  =-!"// (10X,"F:  "  ,  E 1  5 . 7  , 1  X ,  "L  fi  H  (  +  )  ^  *‘,F15.7/)) 

c  7  r  r-<  m  a  t  (  ”  r?  “  ,  F  1 F  .  7  , "  C.vOLF  CEL  fiY  (-)  =  ",  15.7) 

ot:  rorMATi/  (  1  -«v,  "  r=  ",  F15. 7 , 1C.X  ,"GRCUF  EE  LAY  (♦  )  =  **,E15.7/)) 

1 CC  rp'NfiT (7?  10) 

f  M  ^ 
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F  UNf  MON  rfT  MCA) 
rf'MwrN  o,'rL.Tl*r*G*ALl»AL? 
yr  { c a  .  n.  i’.?)  GO  TO  3 
t r  <0?  .  1  T  .  -  1  7  •  3  )  GO  TC  4 
roTM  =1  ./TfNH(CA) 

r. n  to  ? 

3  “  OT  >-=  1 . 

r.n  T  r  ? 

*  rriTHr-1  . 

p  ■'  ONTTNIjF 

7.-  -f 
CM0 


r'JMOTTOM  T  (OA  ) 

COMMON  n,  ctl.t1  ,S,C,AL1,AL2 
r£-UCfi*PL 

IF  (GAEL  .O'.  17.3)  GO  to  43 

IF  { Cf c L  .IT.  -17.3)  GO  TO  >4 

T  =  (  AL  2+TANMC  ACL  )  )/<  AL1-T  ANH  (C  AEL) » 

GO  tc  u5 

43  T  -  (  AL2+  i.  ) /(  At.  l-i .) 

r,  o  T0  4^ 

h 4  T  =  ( Al  2-1. ) /( AL 1+1. ) 

4  5  OONTINIJF 
tup K 
r  m  o 


FUNCTION  FT  (  c  a  ) 

common  n,  ~l. Ti  ,G,0,ALl,Al_? 

FTrryc(  C  ’  C  T  *  C )  ‘  T  (CA  I"  (AL  1  +  00TH  (CDMl )  H  .  ) 
v  +  Cl.-AL?-  CCTH<CA*Tt)>*t-*XP<-C*CA*C) 

P  F  3  u  F  N 

r  Njn 
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COMPLETED  CASES 


In  this  section  there  are  presented  graphical  results  produced 
by  the  computer  programs  for  various  cases. 

In  figures  2-12  are  the  results  of  the  basic  theory  for  one  set 
of  parameters,  omitting  conduction  loss.  Non-apodized  independent 
conductors  are  considered.  Graphs  are  presented  for  wave  number 
(plus  wave),  group  delay  (plus  wave),  wave  number  (minus  wave), 
group  delay  (minus  wave),  normalized  dispersion  (plus  wave),  radiation 
resistance  (plus  wave) ,  radiation  resistance  (minus  wave) ,  total 
radiation  resistance,  total  radiation  reactance,  negative  of  inser¬ 
tion  loss  (plus  wave)  and  negative  of  insertion  loss  (minus  wave). 

In  figures  13-16  are  the  results  of  the  microstrip  model  for 
the  same  set  of  parameters  as  above.  Presented  are  graphs  of  input 
resistance,  input  reactance,  negative  of  insertion  loss  (plus  wave) 
and  negative  of  insertion  loss  (minus  wave) . 

A  second  set  of  parameters  is  considered  in  Figures  17-28.  The 
basic  theory  is  employed  for  non-apodized  independent  conductors. 

There  are  presented  graphs  for  radiation  resistance  (minus  wave) , 
radiation  resistance  (plus  wave),  total  radiation  resistance,  total 
radiation  reactance,  negative  of  insertion  loss  (minus  wave)  and 
negative  of  insertion  loss  (plus  wave).  Results  are  obtained  for 
N=1 ,  N=2,  N=8  and  N=100. 

Another  set  of  parameters  is  considered  in  all  of  the  figures 
29-46. 

In  figures  29-31  the  radiation  resistance  (minus  wave)  is  pre¬ 
sented  for  the  cases  of  no  apodization,  apodization  in  strip  width 
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and  apodization  in  center  to  center  spacing.  The  basic  theory  for 
independent  conductors  is  here  considered. 

In  figures  32-34  there  are  presented  graphs  for  radiation 
resistance  (minus  wave)  for  the  cases  of  no  apodization,  apodization 
in  strip  width  and  apodization  in  center  to  center  spacing.  Here 
the  basic  theory  for  normal  modes  has  been  considered. 

Figure  35  presents  the  radiation  resistance  (minus  wave)  and 
radiation  resistance  (plus  wave)  for  the  basic  theory  with  indepen¬ 
dent  conductors. 

Figure  36  presents  the  radiation  resistance  (minus  wave)  and 
radiation  resistance  (plus  wave)  as  above,  with  no  ground  planes. 

In  figures  37-38  are  presented  the  radiation  resistance 
(minus  wave)  and  radiation  resistance  (plus  wave)  for  the  basic 
theory  with  normal  modes;  for  the  fundamental  mode  and  for  n=3. 

In  figure  39  the  radiation  resistance  (minus  wave) ,  radiation 
resistance  (plus  wave)  and  total  radiation  resistance  are  presented 
for  the  basic  theory  with  independent  conductors  and  no  ground 
planes,  for  N=l. 

In  figure  40  the  radiation  resistance  (minus  wave)  and  radiation 
resistance  (plus  wave)  for  the  case  as  above,  with  N=3,  are  presented. 

In  figure  41  the  radiation  resistances  (plus  wave)  are  presented 
for  the  basic  theory  with  independent  conductors,  for  N=l,  N=2,  N=3 
and  N=4 . 

In  figure  42  the  radiation  resistances  (plus  wave)  for  the  basic 
theory  with  independent  conductors,  for  N=4,  are  presented  for  three 
different  gap  thicknesses. 


In  figures  43-44  there  are  presented  the  radiation  resistance 
(plus  wave)  and  radiation  resistance  (minus  wave)  for  the  basic 
theory  with  normal  modes,  for  N=32. 
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